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We present a parameter estimation technique based on performing joint measurements of a weak 
interaction away from the weak-value-amplification approximation. Two detectors are used to collect 
full statistics of the correlations between two weakly entangled degrees of freedom. Without dis¬ 
carding of data, the protocol resembles the anomalous amplification of an imaginary-weak-value-like 
response. The amplification is induced in the difference signal of both detectors allowing robustness 
to different sources of technical noise, and offering in addition the advantages of balanced signals for 
precision metrology. All of the Fisher information about the parameter of interest is collected. A 
tunable phase controls the strength of the amplification response. We experimentally demonstrate 
the proposed technique by measuring polarization rotations in a linearly polarized laser pulse. We 
show that in the presence of technical noise the effective sensitivity and precision of a split detector 
is increased when compared to a conventional continuous-wave balanced detection technique. 


Introduction Anomalous amplification [1] has been 
shown to be advantageous for precision metrology. Such 
an amplification provides a way to increase a signal while 
decreasing [2] or retaining the technical-noise floor [3, 4]. 
As a result, the sensitivity and precision of measurements 
limited by technical noise can be effectively improved, fa¬ 
cilitating the saturation of the standard quantum limit. 
Anomalous amplification was first proposed for metrol¬ 
ogy with the introduction of the Weak Value (WV) of 
an observable [1, 5], and parameter estimation proto¬ 
cols defined after it are usually known as Weak-Value- 
Amplification (WVA) techniques. The WV of an ob¬ 
servable is obtained by post-selecting the state of a sys¬ 
tem after a weak interaction with a meter system. In 
WVA, such measurements in the system induce a dis¬ 
carding of data counts in the measurements of the meter. 
In addition to the notion that the state of the system 
is post-selected after the weak interaction, we consider 
post-selection as the process of selecting and process¬ 
ing desired events, which, for WVA, results in discarding 
data in the meter. Due to the interference of the pre- and 
post-selection states of the system the WV can take large 
complex values outside the eigenvalue spectrum of the 
observable, which defines the anomalous amplification in 
WVA. Discussion about the quantum interpretation of 
such a phenomenon can be found in Refs. [6-9]. Many 
recent applications of WVA for metrology have been done 
in classical optics, where the interference can be under¬ 
stood using standard wave mechanics [10, 11]. 

Strong postselection is necessary for anomalous ampli¬ 
fication in WVA techniques, but discarding data counts 
has been the target of criticism, and even considered 
“harmful” for metrology [12-14]. However, it has been 
shown, theoretically and experimentally, that the statis¬ 
tical information collected by the measurements is in¬ 


significantly reduced because the amplified signal can 
compensate for the reduced detection flux resulting from 
postselection [3, 4, 15-19]. Such a result is possible un¬ 
der an almost orthogonal pre- and post-selection proce¬ 
dure in the system allowing one to collect nearly all of 
the Fisher information using only a small subensemble 
of measurement counts in the meter [3]. For example, 
the shot-noise limit defined by the input number of pho¬ 
tons used in measurements of the small velocity of one of 
the mirrors of a Michelson interferometer can be reached 
using WVA [16]. Moreover, it was recently shown that 
by postselecting only 1% of the photons, when measur¬ 
ing small optical deflections constrained to intrinsic elec¬ 
tronic detector noise, 99% of the total available Fisher 
information can be recovered [4]. 

Signal amplification while avoiding detector saturation 
sparked interest in WVA as a precision metrological tech¬ 
nique several years ago [20-22], and was recently shown 
to be essential for the technical-noise mitigation advan¬ 
tages [3, 4]. We will introduce the possibility of inducing 
anomalous amplification without the need of discarding 
data and without any loss of Fisher information. 

Striibi and Bruder [23] recently proposed a precision 
measurement technique for measuring time delays of light 
by carrying out a full measurement of the two weakly cor¬ 
related degrees of freedom: frequency and polarization. 
They concluded that even for low-resolution detectors, 
the scheme is robust against systematic errors and fluctu¬ 
ations in alignments of the experimental setup. We show 
that by subtracting the readouts of the two detectors in 
such a measurement, a WVA-like response is obtained 
in the difference signal. This anomalous amplification 
behaviour clarifies and extends the results reported in 
Ref. [23], since besides having the amplification response 
typical of the WVA approach (without discarding data), 
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our protocol adds the benefits of balanced detection for 
precision metrology. We show that this technique allows 
us to recover all of the Fisher information of the esti¬ 
mated parameter. In addition, it permits the removal of 
systematic error in the measurement of the shift in the 
difference signal by tracking the unshifted sum signal as 
well. 

Theoretical Framework Following the formalism as in 
WVA, we describe the unitary evolution of a two-party 
system as 1/ = exp {—igq^A), where A is a binary de¬ 
gree of freedom (qubit) controlling the encoding of the 
information about the interaction parameter g in the con¬ 
tinuous (meter) observable q [24]. In contrast to WVA, 
where a small set of measurements for q are taken into 
account due to postselection, the operator q after the 
interaction is always measured and conditioned to one 
of two detectors. The proposed procedure is done by 
preparing the initial global state as the product state 
kin) ® with = (|0) + |l))/\/2, and by track- 

ing q conditioned to the measurement basis |Ti^ 2 ) = 
(|0) ± ie'^^\l))/^/2 on the qubit/system. Here |0) and 
|1) are the eigenvectors of A, where A = |1)(1| — |0)(0|, 
and projections to the two detectors are labelled by 1 
and 2. The small phase e defines the measuring basis 
on the equatorial plane of the Bloch sphere, where the 
initial prepared state for the qubit also lies. The 

probability distributions measured on the detectors take 
the form 

= ^[l^sm{e + 2gq)]P{q), (1) 

where P{q) = \i^{q)f. We will use a Gaussian state ip{q) 
with variance for the preparation in As an example. 
Fig. 1 shows a schematic of the measuring technique for 
an optical setup, where the observable A is represented 
as the which-path degree of freedom in an interferometer, 
and the phase e controls the interference. 

Under the assumption of a weak interaction, i.e. 
2ga <C min {1, tan e}, we can express 

sin(e + 2gq)P{q) ^ sin(e) P (q — 2ga^ cot e) . (2) 

The peak value of the distribution of Eq. (2) is smaller 
by a factor sine and the position of the peak is shifted 
by an amount Sq = 2ga^ cot e with respect to P{q). The 
sum and difference of the distributions take the form [25] 

P+(g) = Pi+P2 = P(g), (3) 

P- {q; g) = P 2 - Pi ^ Sin(e) P{q- 2ga‘^ cot e) . (4) 

After a large number N of independent measurements of 
g, q^^^ = gi,...,gAri on detector 1 and q^‘^^ = gi,...,gAr 2 
on detector 2, the sum distribution (3) reproduces the 
quantum probability distribution of the input state, and 
the difference distribution (4) has a shifted attenuated 
peak similar to WVA. Note that the measured shift is 



FIG. 1. A/mosTbalanced weak values detection technique for 
measuring a small parameter g in an optical setup. The states 
|0) and |1) correspond to the two paths in the interferome¬ 
ter, a small unbalacing phase e controls the interference, and 
estimation of g is obtained from the sum and difference of 
the two measured distributions. The red solid lines represent 
the directly measured distributions Pi and P 2 , and the blue 
dashed lines are the distribution components as in Eq. (1). 


in the difference probability distribution and not in the 
wave-function '0(g) itself, as it is the case in WVA. In 
fact, the weak values for the measurements are given by 
i, 2\A\^in)/{^1,2\^in) = COSe/(1 T sine) ~ 

=F^(1 ± e) for e ^ 1, and no (anomalous) large weak 
value is induced. Estimations of averaged values for e 
and g under the weak interaction approximation can be 
obtained as 


e = sin 



( N^-NA 
VV1+V2; ’ 
— (g)+] tane 


( 5 ) 

(6) 


where {q)± = (E*^"i qf^ ± ^) / (A ^2 ± A^i), and 

is the measured variance of P+(g). By making e ^ 1 
and preparing a large variance input state, a large shift is 
induced and small values of g resolved. This behaviour is 
similar to the amplification of the WVA technique with 
an imaginary weak value [6], and the quadratic response 
of the WVA post selection probability with respect to e, 
sin^(e/2) ^ e^/4, is replaced with a (larger) linear re¬ 
sponse of the difference signal, sine ^ e in Eq. (4). 

The maximum amount of information about the pa¬ 
rameter g that can be extracted from the measurements 
is given by adding the Eisher information for both detec¬ 
tors [26], 


Fg = Fi + F 2 = (7) 

where Fj = Ni ^ (^^) 

cient estimator for g, which saturates the Eisher infor¬ 
mation (7) in absence of noise. The inverse of the Eisher 
information is known as the Gramer-Rao bound or GRB. 
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This bound is the smallest possible variance of an unbi¬ 
ased estimation of the parameter g. The smallest pos¬ 
sible standard deviation for measurements of g is given 
by Ag^^^ = Fg = 1/ (^2^/Na^^ and any source of 
noise would increase it. This result shows the standard 
quantum limit or shot-noise dependence with respect to 
V^, characteristic of N independent measurements. 

Comparison to other techniques We now consider two 
alternative approaches to measure the parameter g^ and 
compare them to our proposed protocol. We note that all 
three techniques are upper bounded by the same Fisher 
information, . 

The first approach is formulated by noticing that 
the unitary evolution corresponds to a translation of 
g in the canonical momentum of g, and that the an¬ 
cillary system A is not necessary. The protocol con¬ 
sists on preparing N copies of '0(g), applying the evo¬ 
lution Ust = exp(—igg), and performing measurements 
of the conjugate momentum p (instead of g). The vis¬ 
ibility of the shift, i.e. the ratio of the shift and the 
standard deviation of the measured state, is given by 
6plap = gl{ll2a) = 2ga. This approach is known as the 
“standard technique” when similar comparisons to WVA 
are introduced. The visibility of our technique takes the 
form 6q/a = {2ga‘^ cot e)/a = 2ga cot e, giving an advan¬ 
tage of ^ 1/e for e ^ 1 when measuring small values of g. 
This amplification is an important key on how the WVA 
and our protocol are superior in technical-noise-limited 
experiments. 

The second approach is the conventional WVA tech¬ 
nique, where measuring g, using only one detector, is con¬ 
ditioned to the postselection |Tj) = (|0) — e*^|l)) /^/2. 
This procedure is equivalent to removing the balanc¬ 
ing phase 7t/2 in Fig. 1 and tracking only the dark 
port of the interferometer. The probability distribu¬ 
tion for such a measurement is given by P^^^(g) = 
sin^ (e/2)P(g —2^cr^cot(e/2)). So, even though the shift 
of the peak of our protocol is half the shift of the peak 
in the WVA technique for a given e (for e ^ 1), the peak 
value of P-(g) is much larger than in P^^^(g) since 
e ^ e^/4. This result plus the background noise sub¬ 
traction characteristic of differencing signals allow us to 
experimentally induce smaller possible values of e than in 
WVA, which offers technically advantageous larger am¬ 
plification than in WVA. 

The Fisher information for the WVA technique, un¬ 
der the weak interaction approximation, is given by 
pWYA _ 47Vcr^ cos^(e/2) [26]. WVA measurements 
asymptotically recover all of the Fisher information, given 
by 4Vcr^, only in the anomalous weak value regime 
{2ga ^ e/2 ^ 1) [3], but the proposed differencing 
technique collects always all of the information. WVA 
and the almost-balanced weak values technique offer sim¬ 
ilar amplification behaviour, and noise mitigation advan¬ 
tages will depend strongly on the metrological specific 


Piezo-driven HWP 

(45° (f) (jJot) 



FIG. 2. Experimental setup for measuring small angular ve¬ 
locities ujQ of a piezo-driven half wave plate (HWP). The angu¬ 
lar rotation induces changes of polarization in the laser field. 

experimental task. For example, WVA could still be a 
favourable technique in situations where detector satura¬ 
tion is the predominant limiting factor. 

Finally, note that both of the above techniques (stan¬ 
dard and WVA) require prior measurements of the un¬ 
shifted peak and variance of the input distribution P(g). 
Ours does not, since the sum distribution P+(g) offers 
these measurements simultaneously. This is one of the 
advantages of using two detectors instead of just one. 

Experimental implementation As a proof-of-principle 
of the technique in the optical domain, we performed 
measurements of small polarization changes in a linearly 
polarized pulse of laser light. A piezo-driven half-wave 
plate (HWP) played the role of the interferometer in 
Fig. 1, where we used polarization instead of the which- 
path degree of freedom, and time as the variable g (with 
a = r). The piezo actuator rotated the HWP in time 
by an angle 0 + cjot, which defined the tunable phase as 
e = 40, and the angular velocity of the rotating HWP as 
the parameter of interest, g = 2uJo (see Fig. 2 and supple¬ 
mental material for experimental details). The estimates 
for ujQ using a split detector showed an almost perfect lin¬ 
ear response and small standard deviations even without 
the need of pulse averaging. The angle 0 was 4.972(5) 
mrad on average. This unbalanced phase is significantly 
smaller than previously reported postselection angles us¬ 
ing WVA techniques, proving the possibility of larger am¬ 
plification under our proposed protocol. 

Besides the two (standard and WVA) experimental 
techniques mentioned before, we compare here our pro¬ 
tocol to a conventional balanced optical technique for 
measurements of cJo, where no ancillary system nor mea¬ 
surements of peak shifts in a distribution are required. 
By replacing the Gaussian pulse with a continuous wave 
(c.w.) beam, the perfectly balanced (0 = 0) signal takes 
the form /_ (t) = Iq sin(u;o^) and the value ujq can be 
recovered. The measured intensities are very small com¬ 
pared to our proposed technique, since ujot 0. Thus, 
the technical advantages of our technique rely principally 
on time shift measurements of a Gaussian profile with 
controllable amplitude (/osin40), instead of measuring 
very small voltage amplitudes as it is the case of the c.w. 
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conventional balanced technique. For example, in order 
to obtain a visible signal (signal-to-noise ratio slightly 
larger than one) using a c.w. beam, 1024 periods of in¬ 
tegration time were required. Our technique gives a bet¬ 
ter signal-to-noise ratio (> 29) even without the need of 
pulse averaging [26]. 

The best possible variance of our measurements are 
defined by the CRB for cjq- Following equation (7), 


A, jCRB _ ^ ^ 

4v^r’ 


(8) 


where a is the number of averaged pulses. The exper¬ 
imental standard deviations of the measurements with 
the split detector were estimated to be between 20 and 
37 times the calculated ones using Eq. (8). Such devia¬ 
tion is small considering that no frequency filters, lock- 
in amplifiers or any other electronic processing device 
was used. In order to have good photon number statis¬ 
tics and a clearer raw signal, the split detector of Fig. 2 
was replaced with two single photon counting modules 
(SPCM). The input power was attenuated before the 
Acoustic Optic Modulator (AOM) so that the peak de¬ 
tected photon rate was slightly smaller than 10^ counts/s 
on each SPCM. Fig. 3 shows the standard deviation on 
the estimation of ujq as a function of the number of aver¬ 
aged pulses a, and its comparison to eq. (8). The system 
follows the expected shot-noise behaviour, and 

it has (inset of Fig. 3) a precision close to the shot-noise 
limit. There is a slight decrease in the relative precision 
(Aujq/Aujq^^) with the increase on the number of aver¬ 
aged pulses, which is due to long-term drift [26]. Fig. 3 
shows results for a total collection time of up to 5000 
sec (50x100 pulses) and a factor of no larger than ^1.6 
away from the shot-noise limit, in comparison to a factor 
of ^1.2 for shorter acquisition times. 

Conclusions We have introduced a metrological tech¬ 
nique for measuring small interaction parameters based 
on almost balancing both signals of a joint measure¬ 
ment. Even though the protocol does not require the 
postselection of WVA techniques, anomalous amplifica¬ 
tion is recovered. A small phase controlling the unbal¬ 
ancing of the difference signal plays a similar role to the 
postselection angle in WVA. By subtraction of the sig¬ 
nals, background noise floor is removed and also a lin¬ 
ear response with respect to the small phase e is ob¬ 
tained. Such linear response is larger than the quadratic 
response of WVA techniques, and it is the principal ad¬ 
vantage of our technique by allowing larger amplification 
(smaller e) in technical-noise-limited experiments than 
in WVA techniques. In addition, no prior knowledge of 
the peak and the variance of the distribution of the in¬ 
put state is required, so systematic error-free estimation 
of the parameter of interest is possible. We believe our 
technique will find interesting applications on precision 
metrology, within and outside optical systems. The use 



FIG. 3. Standard deviation on estimations of ojq as a function 
of the number of pulses averaged, a. Pulses with r = 0.1 s 
and N = 5.246(2) x 10^ were used. The driving signal on 
the piezo-mounted HWP was 60 V peak-to-peak at A = 1 
Hz. Each point (blue) was obtained from the statistics of 
100 measurements, and the solid line (red) is given by the 
shot-noise (eq. 8). The inset shows the ratio between the 
experimental standard deviation and the shot-noise. 


of post-selection required for WVA has been challeng¬ 
ing in experimental set-ups outside of optics. Our pro¬ 
posed almost-balanced weak values technique overcomes 
the necessity of post-selection and might have direct ap¬ 
plications in non-optical experiments. 

Anomalous amplification in our differencing technique 
allows us to collect all of the Eisher information that 
WVA techniques collect only in the optimal asymptotic 
limit of a large/anomalous weak value. This result proves 
that anomalous amplification can be obtained from the 
weak system-meter interaction assumption alone, where 
strong measurements in the system condition the reading 
in the meter. Technical advantages of our proposed tech¬ 
nique and WVA over standard techniques depend upon 
such a coupling to an ancillary system. 

A recent related work reported experimental advan¬ 
tages of the original proposal of Ref. [23] with re¬ 
spect to the WVA protocol when estimating ultra-small 
phases [27]. 
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SUPPLEMENTAL MATERIAL 


CALCULATION OF THE FISHER INFORMATION 


Fisher information is the maximum amount of information that can be extracted about an unknown parameter {g 
in our case) from a random variable q. It is calculated as the variance of the score. The score is a measure of the 
sensitivity of the probability distribution P{q; g) with respect to the parameter g^ and it is given by 


V{q-,g) = ^^\nP{q-g). 

In our case, we have two different probability distributions, coming from two detectors, 

1 




[lTsin(e + 25g)] , 


2\/^a 

where i = {1,2}. The mean value of events collected by each detector is given by 


N / 2 2 \ 

Pii9;q)dq=- (l T - sine) 


(9) 


( 10 ) 


( 11 ) 


and A = + A 2 is the total number of measurements, which (in a lossless system) is equal to the number of prepared 

events. 

Note that the mean of the score vanishes, 


(V) = I Pi{q;9) 

Pi{q-,g)dq 


d 

— \nPi{q-,g) 


dq 


y 
dg \_. 


= 0 , 


( 12 ) 
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so the final expression for the Fisher information of g takes the form 

Fg = F^+F2= N^iy^) + N2{Vi) 

= 4Ny, 


(13) 


where 


Fi = Ni 


L 


-oo-P*(9;5)L dg 


dPi{q;9) 


dq = 2Ny 1 ± (1 - iff 






(14) 


The result in eq. (13) is general for any value of g and e, meaning that it is always possible to recover all of the Fisher 
information by using both detectors. The shift in the difference signal is obtained only under the weak interaction 
approximation, 2ga <C min {1, tan e}. To induce the anomalous amplification is also required the almost-balanced 
condition, e ^ 1. In this sense, our proposed technique lies in the regime of 2ga e 1. 

Now we compare the result of Eq.( 13) to the Fisher information when using the WVA technique. The intensity in 
the dark port under the weak interaction approach is given by 

= y- sin" (I + m) (15) 

V27rcr / 


and after a straightforward calculation we get the number of detected events and the Fisher information, 

N r. 


WVA 


N 


pWVA ^ 2N(j‘^ 




cose 


1 + (1 — 4^^cr^)e ^ cose 


Nsm\e/i) + 0[gf, 

mA cos^{e/2)+ 0[g\^, 


respectively. 


(16) 

(17) 


EXPERIMENTAL METHODS 

An external cavity 795 nm wavelength diode laser and an Acoustic Optical Modulator (AOM) are used to prepare 
pulses with an intensity profile I{t) = Jq exp (—t^/2T^). A Gian-Taylor calcite polarizer prepares the polarization 
state in \D) = {\H) + jE)) /a/ 2, and a piezo-driven half wave plate (HWP) rotating in time by an angle (/)-{-ujot away 
from the \D) axis changes the polarization of the pulse. We are interested in estimating the angular velocity ujq of 
the rotating HWP (see Fig. 2 in the main document). 

A Wollaston prism was used to separate the horizontal {\H)) and the vertical (|H)) polarization components, which 
were independently measured using a segmented quadrant detector. The sum and difference signals take the form 

(18) 

I-{t) = Jo sin(4(/)) e-(*-‘^o^V 0 )V 2 T^^ 

assuming that coqt ^0^1. These pulses were recorded and the angle <p and the time shift 6t = obtained 

from numerically fitting the distributions. We estimated <p and c^o by using eqs. (5) and (6) of the main document as 
well, nevertheless both estimators gave us identical results within the experimental errors. 

A triangle function in voltage with frequency /^, peak-to-peak voltage Vpp and 60% duty cycle as an input for 
the piezo actuator induced a polarization rotation ramp given by c^o = (10/6)QfV/p/r, where a is the HWP-mounted 
piezo response. Gaussian pulses with a peak power of ^ 136 /iW were sent in phase with the ramp. Fig. 4 shows the 
measurement results of ujq for three different values of fr as a function of Vpp. The technique’s estimates show an 
almost perfect linear response and small error bars even without the need of pulse averaging. The angle <p was 4.972(5) 
mrad on average. The values for the pulse length r obtained from numerically fitting the measured sum intensities are 
5.377(2) ms for 20 Hz, 10.695(3) ms for 10 Hz, and 21.398(7) ms for 5 Hz. In Fig. 4 we use the errors propagated from 
the time shift measurements (error of the mean) which were obtained from the statistics of 60 measured time shifts 
for each point. The experimental standard deviations of the measurements of Fig. 4 were estimated to be between 20 
and 37 times the shot-noise (calculated from Eq. (8) in the main document). 

We were limited by the resolution of the HWP. We measured ^30 /irad (Vpp = 15 V) as the smallest reliable total 
rotation induced by the HWP on the pulses, so we boosted the effective sensitivity of our setup for measuring ujq 
by decreasing the driving frequency on the piezo, fr. Keeping the peak-to-peak voltage on the piezo at 15 V we 
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FIG. 4. Measurements of ujq as a function of the peak-to-peak voltage applied to the piezo actuator. The estimated averaged 
number of collected photons are kept nearly constant, ~ 2.9 x 10^^, by measuring one single pulse when fr =5 Hz, averaging 
two pulses when fr =10 Hz, and averaging four pulses when fr =20 Hz. The averaged Signal-to-Noise ratios are between 29 
for Vpp = 15 V up to 190 for Vpp = 75 V. 


set fr = 2.5 mHz, r = 30.08(1) s and (f = 4.976(6) mrad, and measured our lowest reported angular velocity of 
ujQ = 161 ± 22 nrad/s. Our reported net rotation of 32 ± 4/irad in a rotation time of 200 seconds was obtained by 
averaging sets of five pulses and having 20 statistical points. Our measurements for these long periods were affected 
by long-term drifts due to temperature and air pressure fluctuations, so 100 consecutive pulses (out of 360 collected) 
were selected such that the angle 0 , the peak laser power and the measured time shift remained relatively constant. 


POLARIZATION/PHASE ROTATIONS 


We prepare the system in the state 

|i>o) = l^-o) 0 E{t) = T (\H) + \V)) 0 Eo 


( 20 ) 


where Eq is the peak electric field of the Gaussian pulse. We use a classical description for the long non-Fourier-limited 
pulses of light used in our experiment, but the theory can also be applied to coherent pulses of light or single photons. 
After passing through the piezo-driven HWP, the polarization is rotated as 


l^o) ^ —j=^ [{cos(20 + — sin(20 -h 2ujQt)} \H) 

v2 

T |cos(20 T T sin(20 T 2c(Jo^)} |f^)] ^ ^ ^ 


( 21 ) 


and after separating the two polarization components, using a Wollaston prism, the recorded intensities in the split 
detector are given by. 


/i = -^ [cos(2(/> + 2 i 0 ot) — sin(2^ + 2u}ot)]^ e~* 

(22) 

= h [1 _ sin (4(^ -(- 4wot)] e“* , and 

(23) 

/2 = h [1 + sin (4<)) + 4woi)] . 

(24) 

Finally, calculating the sum and difference signals we get 

/+ =/i +/2 = lo and 

(25) 

/_ = /2 - /i = lo sin {4:<p + Acoot) 

(26) 

«/osin(40) 

(27) 

« /osin(40) 

(28) 


where ujqt <C 0 <C 1 . 
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Note that if we change the basis in equation (20) from \H) and \V) to circularly polarized states, 

\R) = ^i\H)-i\V)),e.nd 
\L) = ^{\H) + i\V)), 

the initial state takes the following representation (up to a global phase) 


\^o) = ^{\R)-i\L))^Eoe-^^/^^\ (31) 

This representation defines an exact mapping of the theoretical protocol presented in Fig. 1 of the main document, 
where both states {\R) and \L)) are separated into the two different paths of the interferometer, and there is a relative 
phase of 4(0 + ujQt) between both arms. In fact, the state after the light pulse passes trough the piezo-driven HWP 
takes the form (eq. [21]), 


gi2(0+u;ot) 1^^ _ ^ g-i2(0+u;ot) 


Under this representation, the left and right circularly polarized components of light split into the two arms of the 
interferometer. If Fourier-limited pulses were used, the quantity 4uJo could be understood as a relative laser frequency 
shift between the two paths. 



